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ABSTRACT 

The main goal of this research is to get better insights into the properties of the 
plasma filled magnetospheres of black holes by means of direct numerical simulations 
and, ultimately, to resolve the controversy surrounding the Blandford-Znajek mech- 
anism. Driven by the need to write the equations of black hole electrodynamics in 
the form convenient for numerical applications, we constructed a new system of 3+1 
equations which not only has more traditional form than now classical 3+1 system 
of Thorne and Macdonald but which is also more general. To deal with the mag- 
netospheric current sheets, we also developed a simple model of radiative resistivity 
based on the inverse Compton scattering of background photons. The results of nu- 
merical simulations combined with simple analytical arguments allow us to make a 
number of important conclusions on the nature of the Blandford-Znajek mechanism. 
We show that, just like in the Penrose mechanism and in the MHD models of Punsly 
and Coroniti, the key role in this mechanism is played by the black hole ergosphere. 
The poloidal currents are driven by the gravitationally induced electric field which 
cannot be screened within the ergosphere by any static distribution of the electric 
charge of locally created pair plasma. Contrary to what is expected in the Membrane 
paradigm, the energy and angular momentum are extracted not only along the mag- 
netic field lines penetrating the event horizon but along all field lines penetrating the 
ergosphere. In dipolar magnetic configurations symmetric relative to the equatorial 
plane the force-free approximation breaks down within the ergosphere where a strong 
current sheet develops along the equatorial plane. This current sheet supplies energy 
and angular momentum at infinity to the surrounding force-free magnetosphere. The 
Blandford-Znajek monopole solution is found to be asymptotically stable and causal. 
The so-called horizon boundary condition of Znajek is shown to be a regularity con- 
dition at fast critical surface. 

Key words: black hole physics: general - magnetic fields - methods: numerical 



1 INTRODUCTION 

As the result of great advances in observational astrophysics 
during the previous several decades black holes are no longer 
regarded as peculiar solutions allowed by the Einstein equa- 
tions which may or may not have anything to do with real as- 
tronomical phenomena. It is now widely believed that black 
holes are common in the Universe and play a key role in 
the most violent space events such as activity of galactic 
nuclei and gamma-ray bursts. Enormous amounts of energy 
released during such events can have two different origins. 
First of all this can be the gravitational energy of matter re- 
leased during accretion onto an already existing black hole 
or during the gravitational collapse leading to formation of a 
new black hole. On the other hand, this can be the rotational 
energy of a black hole itself. 

The notion of rotational energy of a black hole emerged 
as the result of theoretical discovery made by Penrose(1969) 



who found that a fraction of the Kerr hole mass can be 
converted, at least in principle, into the energy of surround- 
ing matter or radiation (see also Cristodolou 1970). He has 
shown that two particles may interact in the neighbourhood 
of the black hole in such a way that one of the particles 
acquires negative energy and eventually disappears into the 
hole whereas the other particle carries away excessive pos- 
itive energy. Unfortunately, as this was shown later, such 
interactions are rather special and must be very rare under 
typical astrophysical conditions, rendering the Penrose pro- 
cess inefficient (Bardeen et al. 1973). However, it was also 
found that the electromagnetic field can be used to extract 
the rotational energy of black holes too. 

First, Goldreich and Julian(1969) analysed the vacuum 
solution for a rotating neutron star with a dipolar magnetic 
field aligned with the rotational axis. They argued that the 
rotationally induced electric field was strong enough to pull 
charged particles from the stellar surface and, thus, fill the 
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surrounding space with plasma. Using the force-free approxi- 
mation to describe the produced magnetosphere they argued 
that an electromagnetically driven wind would carry away 
rotational energy and angular momentum of the star. 

Then, Wald(1974) found a rather interesting particu- 
lar solution of the vacuum Maxwell equations in the Kerr 
spacetime. Far away from the hole this solution described 
a uniform magnetic field aligned with the rotational axis of 
a black hole. However, near the black hole it described a 
strong electric field as well. Moreover, just like in the prob- 
lem of an aligned rotator this "gravitationally induced" elec- 
tric field had a significant component along the magnetic 
field. Bisnovatyi-Kogan and Ruzmaikin (1976) argued that, 
in the case of astrophysical black holes, a rather strong mag- 
netic field of such kind can be generated in their accretion 
discs. 

Finally, Blandford & Znajek (1977) realised that the 
similarity between the vacuum solution for a Kerr black hole 
and the vacuum solution for a rotating neutron star meant 
the possibility of electromagnetically driven wind from a ro- 
tating black hole, provided the space around the black hole 
could be filled with plasma. Moreover, they argued that, 
under the typical astrophysical conditions, the vacuum so- 
lutions were, in fact, unstable to cascade pair production, 
ensuring a plentiful supply of charged particles. Then they 
developed a general theory of force-free steady-state axisym- 
metric magnetospheres of black holes and found a perturba- 
tive solution for a slowly rotating black hole with monopole 
magnetic field. The key element of this solution was Zna- 
jek's "boundary condition" (Znajek 1977) imposed on the 
event horizon. As expected, this solution exhibited outgo- 
ing electromagnetic fluxes of energy and angular momen- 
tum. Moreover, the electromagnetic mechanism seemed to 
be very robust and the estimated power of the wind was 
high enough to explain the energetics of radio galaxies and 
quasars. 

In their analysis, Blandford & Znajek (1977) used co- 
variant equations of electrodynamics and operated with 
components of the electromagnetic field tensor and four- 
potential. Later, Thorne and Macdonald (Thorne & Mac- 
donald 1982; Macdonald & Thorne 1982) developed a 3+1 
approach, where the equations of black hole electrodynam- 
ics were written in more or less traditional form in terms 
of the spatial vectors of electric and magnetic field as mea- 
sured by the so-called local fiducial observer (FIDO). In this 
work they adopted the system of coordinates due to Boyer 
and Lindquist(1967) which has a number of useful proper- 
ties but also one important drawback. Just like the system 
of Schwarzschild coordinates, it is singular at the event hori- 
zon. As the result, the event horizon appears as a peculiar 
inner boundary of of physical space which required rather 
special treatment. 

A number of authors studied the properties of electro- 
magnetic field near the event horizon (Hanni & Ruffini 1973; 
Hajicek 1974; Znajek 1977,1978; Damour 1978) and gradu- 
ally a picture emerged, according to which the event hori- 
zon could be treated as a rotating conducting surface with 
surface charges, surface currents, and a finite surface resis- 
tivity. This perfectly suited the quest of Thorne and Mack- 
donald for a new formulation of the black hole electrody- 
namics which would make it look similar to the classical 
electrodynamics. Surprisingly enough, the drawback of the 



Boyer-Lindquist coordinates seemed to turn into an advan- 
tage. This theory of the event horizon have made a great 
impact on the current perception of the Blandford-Znajek 
mechanism which is now widely associated with a mental 
picture of magnetic field lines originating from the horizon 
and being torqued by its rotation. It stimulated the develop- 
ment of The Membrane Paradigm (Thorne et al.1986) where 
the event horizon is attributed with a whole range of physi- 
cal properties. However, one has to admit that, in spite of all 
its attractive simplicity and mathematical correctness, this 
construction is purely artificial. Moreover, it is rather wor- 
rying that the emphasis put on the role of the event horizon 
makes the Blandford-Znajek mechanism completely alien to 
the Penrose mechanism, where the key role is played by the 
black hole ergosphere. 

The electrodynamic mechanism together with the hori- 
zon theory is now widely accepted by the astrophysical com- 
munity. In great contrast to this mainstream trend, Punsly 
and Coroniti (Punsly & Coroniti 1990a; Punsly & Coro- 
niti 1990b) and later Punsly (see the review in Punsly,2001) 
completely rejected both these theories. They argued that 
the event horizon cannot be regarded as a unipolar inductor 
because it is causally disconnected from the outgoing wind. 
Indeed, both the fast and the Alfven waves generated at the 
event horizon can propagate only inwards and cannot effect 
the events in the outer space. The apparent lack of a proper 
unipolar inductor in the Blandford-Znajek solution and its 
reliance on Znajek's boundary condition made Punsly and 
Coroniti to conclude that this solution is nonphysical and 
structurally unstable. They developed completely different 
MHD models which seemed to be based on clearer physical 
ideas. In brief, they argue that gravity forces magnetospheric 
plasma to rotate inside the black hole ergosphere in the same 
sense as the black hole and that the magnetic field exhibits 
a similar rotation because it is "frozen" into this plasma. 

Koide (2003) carried out MHD simulations of dipolar 
black hole magnetospheres with plasma energy- density only 
1 + 2 order of magnitude smaller than the energy density of 
electromagnetic field. Although the numerical solution did 
not reached a steady state in these simulations, the obtained 
results seemed to indicate that in this regime the inertial ef- 
fects accounted for approximately half of the extracted en- 
ergy. However, in the black hole magnetospheres filled with 
plasma via pair production the characteristic number den- 
sity of particles is given by the so-called Goldreich- Julian 
density (Goldreich & Julian 1969; Beskin et al.1991; Hirotani 
& Okamoto 1998). Under the typical conditions of super- 
massive black holes in AGNs, this corresponds to the mass 
density of pair plasma hardly exceeding 10~ 13 of the energy 
density of the electromagnetic field. Thus, even if the par- 
ticle density is several orders of magnitude larger than the 
Goldreich- Julian one, the inertial effects are still expected to 
be negligibly small. Given the fact that, in such a degenerate 
regime, the full system of relativistic MHD presents a real 
challenge for numerical methods (see Komissarov, 2001a, 
Gammie et al., 2003), it is difficult not to conclude that 
the electrodynamic approach is more suitable. This is one 
of the reasons for the renewed interest towards the theory 
of force-free electrodynamics (Komissarov 2002a; Blandford 
2002). (This does not mean that the pair production is the 
only possible way of plasma supply in the neighbourhood 
of astrophysical black holes and that the inertial effects are 
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always that small. For example, the coronas of black hole ac- 
cretion discs are likely to be filled with dense plasma pulled 
from the disc surface.) 

The recent numerical studies of the force-free magneto- 
spheres of black holes added to the controversy surround- 
ing the Blandford-Znajek mechanism. On one hand, they 
revealed the asymptotic stability of the Blandford-Znajek 
monopole solution (Komissarov 2001b) and, thus, raised 
doubts about the validity of the causality arguments by 
Punsly and Coroniti. On the other hand, the results for the 
dipolar magnetospheres questioned the virtues of the hori- 
zon theory as well. Indeen, they clearly indicated that the 
key role in the electrodynamic mechanism is played not by 
the black hole event horizon but by its ergosphere (Komis- 
sarov 2002b). Moreover, these simulations also revealed a 
certain deficiency of the force-free approximation as a dis- 
sipative current sheet was formed in the equatorial plane 
within the ergosphere. 

This paper is an attempt to understand the nature of 
the electrodynamic mechanism and to resolve the contro- 
versy surrounding it. In Sec. II we show that the 3+1 equa- 
tions of the black hole electrodynamics can be written in a 
more general and somewhat simpler form than in (Macdon- 
ald & Thorne 1982; Thorne & Macdonald 1982). In Sec.III 
the basic results for the force-free magnetospheres are re- 
derived using the new 3+1 formulation. In order to handle 
the ergospheric current sheets, one has to go beyond the 
force-free approximation and consider the resistive electro- 
dynamics. A model of resistivity based on the inverse Comp- 
ton scattering of the background photons is described in 
Sec. IV. In Sec.V we present the results of numerical studies 
of the black hole magnetospheres obtained within the frame- 
work of resistive electrodynamics. The implications of these 
results for the perception of the Blandford-Znajek mecha- 
nism are discussed in Sec. VI. The more technical informa- 
tion about the properties of light surfaces, the properties of 
the Kerr-Schild coordinate system, and the details of our 
numerical method is presented in the Appendix. 

Throughout this paper we adopt ( h ++) signature 

for the spacetime and assume that the Greek indices range 
from to 3 and refer to spacetime tensors whereas the Latin 
ones range from 1 to 3 and refer to purely spatial tensors. 
In addition, we adopt such units that the speed of light and 
4-k do not appear in the equations of electrodynamics. 



2 3+1 ELECTRODYNAMICS OF BLACK 
HOLE 

Following Macdonald & Thorne (1982) we adopt the folia- 
tion approach to the 3+1 splitting of spacetime in which the 
time coordinate t parametrises a suitable filling of spacetime 
with space-like hypersurfaces described by the 3-dimensional 
metric tensor 7^ . These hypersurfaces may be regarded as 
the "absolute space" at different instances of time t. Below 
we describe a number of useful results for further references. 
If {a: 1 } are the spatial coordinates of the absolute space then 



The 4-velocity of the local fiducial observer, 'FIDO', which 
can be described as being at rest in the absolute space, is 



(/3 2 - a 2 )dt 2 + 2/3idx'dt + ■y ij dx i dx j 



(1) 

where a is called the "lapse function" and (3 is the "shift 
vector" . In this paper we consider only such coordinates that 



ra„ = (-a, 0,0,0). 



(3) 



The spatial components of the projection tensor, which is 
used to construct pure spatial tensors, 



7<*/3 = 9a/3 + n a n f3 , 



(4) 



coincide with the components of the spatial metric 7^. Other 
useful results are 



n" = -(1,-/3*), 



9 = n", 

a 



g = -a 7, 
where 



(5) 

(6) 
(7) 



P l = 7 lJ 13 j , g = det , 7 = det 7^ . 

f3* are the components of the velocity of the spatial grid rela- 
tive to the local FIDO as measured using the coordinate time 
t and the spatial basis {di} (Macdonald & Thorne 1982). 

Landau and Lifshitz (1951) proposed a different ap- 
proach to the 3+1 splitting of spacetime which is known 
as the congruence approach. They also showed that within 
this approach the covariant equations of electrodynamics in 
gravitational field can be reduced to a set of 3+1 equations 
very similar to the usual Maxwell equations in matter (e.g. 
Landau and Lifshitz, 1971). 

The 3+1 electrodynamics of black holes within the fo- 
liation approach was developed by Thorne and Macdonald 
(1982) who used as a starting point the results by Ellis(1973) 
on the 3+1 electrodynamics in the congruence language. In 
fact, they derived rather general 3+1 integral equations of 
electrodynamics as well as the differential equations adapted 
to a particular foliation of the Kerr space-time, namely the 
Boyer-Lindquist foliation. This formulation has been used 
successfully in various theoretical studies of black hole elec- 
trodynamics. However, the Boyer-Lindquist foliation has one 
important disadvantage - it leads to the well known coor- 
dinate singularity at the event horizon. This singularity is 
removed in a different less known foliation, namely the Kerr- 
Schild foliation, which otherwise shares many common prop- 
erties with the Boyer-Lindquist foliation (see Appendix B). 
It has been shown that the Kerr-Schild coordinates do not 
only simplify theoretical analysis but also allow to overcome 
a number of problems concerning numerical simulations (e.g. 
Papadopoulos & Font, 1998, Komissarov, 2001). The corre- 
sponding set of 3+1 differential equations of electrodynamics 
is somewhat different. Moreover, in order to construct a Go- 
dunov type numerical scheme one needs a different type of 
integral equations compared to those given in Thorne and 
Macdonald (1982). All these equations can be obtained from 
the 3+1 integral equations of (Macdonald & Thorne 1982). 
However, they can also be derived directly from the covari- 
ant Maxwell equations, and so easily, that this is worth to 
be shown (see also Landau & Lifshitz, 1971; Koide, 2003). 

The covariant Maxwell equations are (e.g. Jack- 
son(1975)): 



dtg a /3 = 0. 



(2) 



V = 0, 



(8) 
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(9) 



where F al3 is the Maxwell tensor of the electromagnetic field, 
*pafj j s tne F arac i a y tensor and I a is the 4-vector of the 
electric current. The most direct way of 3+1 splitting of 
the covariant Maxwell equations is to write them down in 
components and then to introduce such spatial vectors that 
these equations have a particularly simple and familiar form. 
For example, when eq.(8) is written in components it splits 
into two parts: 

• The time part: 

^ft(a^*F") =0, (10) 

• The spatial part: 

-±=a (aV7*i^) + ^=d t (a^*F») = 0, (11) 

If we now introduce the spatial vectors B and E via 

B l = a*F lt (12) 

E l = ^E 3 , E t = ^e I]k 'F' k , (13) 



where e ijfc = ^yeij/c is the Levi-Civita pseudo-tensor of the 
absolute space, then equations (10,11) read 



V B = 0, 

9 t B + Vx£=0, 



(14) 
(15) 



where V is the covariant derivative of the absolute space. 
Notice, that in order to derive eq.(15) we used the station- 
arity condition (2). Similarly, equation (9) splits into 



V D=p, 

-d t D + VxH = J 
where 
D l = aF u , 
FF = j ij Hj, 



H> = ^e ijk F jk . 



p = al , J = al 



(16) 
(17) 

(18) 
(19) 
(20) 



As one can see, these 3+1 equations have exactly the same 
form as the classical Maxwell equations for the electromag- 
netic field in matter. Applying V to eq.(17) and then using 
eq.(16) one obtains the electric charge conservation law 



dtp + V-J = 0. 



(21) 



Similar to any highly ionized plasma, the pair plasma of 
black hole magnetospheres has essentially zero electric and 
magnetic susceptibilities. In such a case, the Faraday tensor 
is simply dual to the Maxwell tensor 



F <*0 
Here 



(22) 
(23) 

(24) 



is the Levi-Civita alternating pseudo-tensor of spacetime 
and eafi^v is the four-dimensional Levi-Civita symbol. This 
allows to obtain the following alternative expressions for 
B, E, D, and H: 



B % = \e ijk F jk 



Et = F it , 
1 

2 f 

Hi = *F U . 



l> = ^e ijk *F jk , 



(25) 
(26) 
(27) 
(28) 



Moreover, from the above definitions one immediately finds 
the following vacuum constitutive equations: 



E = aD + (3xB, 
H = aB - (3xD. 



(29) 
(30) 



At infinity, provided the coordinate system becomes 
Lorentzian there, one has a = 1, (3 = and, hence, 



B = H, 



E = D. 



One can easily obtain the following covariant forms of defi- 
nitions (12,13,18,19,20): 



B» = -*F> 11, n„, 

E" = -\^ v e vaBl k a *F \ 



D" 



p^ 



n v , 



H» = -^e„ aPl k a F fJ \ 

71/ 



(31) 
(32) 
(33) 
(34) 

(35) 
(36) 



where k a = d t - All these spacetime vectors are purely spa- 
tial, that is they have zero time component. 
In terms of four-potential W M one has 



Fnir 



-2Wi 



Introducing the scalar potential $ and the vector potential 
A as 

$ = -U t , (37) 

A'^^Ai, A,^U % (38) 

one obtains the familiar looking results 

E = -V$ - d t A, (39) 

B = VxA. (40) 

The covariant definitions are 

A»='f v U v , $ = -W M fc". 

From the covariant equation of motion for a particle 
with the four- velocity u v , the four-momentum p M , and the 
electric charge q, 

^ = qF,„u v , (41) 
one easily obtains the following energy equation 
%= qE .v, (42) 
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where e p = — pt is known as the "energy at infinity" and 
v" 1 — dx z /dt is the spatial velocity vector. 

Both in the Boyer-Lindquist and the Kerr-Schild coor- 
dinates, where d$g vll = 0, one has an equally simple result 
for the angular momentum at infinity, l p — P4>, 



^=q(E + vxB)-m, 



(43) 



where m — 84,. 

In terms of Maxwell tensor the electromagnetic stress- 
energy-momentum tensor is 



Tt = F^F„ 7 -i(F 75 F^)n. 

In terms of B,E,H, and D we have 



1 

2a 
1 

a 



(E D + B H), 



T\ = --e^ k E j H k , 



(44) 

(45) 
(46) 
(47) 



T\ = -e ijk D j B k 
a 

T\ = ~-(D'E 3 + B'Hj) + -^—(E-D + B-H)8). (48) 
a 2a J 

The covariant energy-momentum equation is 

V V T\ = -F^P. (49) 

Both in the Boyer-Lindquist and in the Kerr-Schild coordi- 
nates, as in any other coordinates with cyclic coordinate <j> 
such that dtg^iv = d^g^v = 0, 3+1 splitting of eq.(49) leads 
to the following energy 

d t e + V-S = -(E-J), (50) 
and the angular momentum 

d t l + V-L = -(pE + JxB)-m, (51) 
equations for the electromagnetic field. Here 

e = -aT\ = ^(E-D + BH) (52) 
is the volume density of energy at infinity, 
I = aT\ = (DxB)-m (53) 
is the volume density of angular momentum at infinity, 
S = ExH (54) 
is the flux of energy at infinity, and 

L = -(E-m)D - (Hm)B + ^(E-D + B H)m (55) 

is the flux of angular momentum at infinity. 

It is well known that evolution equations like (15) and 
(17) can be written as conservation laws. In our case we 
obtain 



d t B l + VjiT 3 = 0, 
and 



where 

K l3 =e ijk E k and V 3 = -e ijk H k . 



(56) 
(57) 
(58) 



are the magnetic and the electric field flux tensors. The cor- 
responding integral equations are now easily obtained using 
the divergence theorem 



K lJ dSj = 0, 



^ / D l dV + [ V 3 dS 3 = - [ TdV, 
Jv J sv Jv 



(59) 



(60) 



1 sv Jv 
where 5V is the closed boundary of the spatial volume V; 
dV = ^/ 7 ydx 1 dx 2 dx i is the infinitesimal metric volume of 
the absolute space, and dSi = eij k dx 3 ^dx k 2 ^ is the infinites- 
imal metric surface element of 8V. One can see that both 
the differential and the integral 3+1 equations of black hole 
electrodynamics are, in fact, very similar to the correspond- 
ing equations of Minkowski spacetime. This suggests that 
many well known techniques can, almost readily, be used to 
solve this equations numerically. 

In (Macdonald & Thorne 1982; Thorne et al.1986) the 
3+1 equations of electrodynamics are written in terms of 
the electric and the magnetic field vectors as measured by 
the local FIDO. We denote this vectors as E and B. From 
definitions (33,31) it follows that 



B = B and E = D. 



(61) 



The main benefit gained by the introduction of FI- 
DOs and, hence, the quantities measured by these observers, 
comes from the principle of equivalence. These observers can 
be treated as locally inertial observers and in their frames all 
local physical phenomena are governed by the laws of spe- 
cial relativity. For example, in order to close the system of 
Maxwell equations we need an additional constitutive equa- 
tion relating the electric current to the electromagnetic field, 
namely the Ohm law. Thanks to the principle of equiva- 
lence, this can be done in terms of quantities measured by 
FIDOs, in the same way as in the special relativistic elec- 
trodynamics. From equation (36) it follows that p is, in fact, 
the electric charge density as measured by FIDOs. As shown 
in (Thorne & Macdonald 1982) the electric current density 
as measured by FIDOs, j, is related to J via 



J = aj — p(3. 



(62) 



The second term in this equation accounts for the motion of 
spatial grid relative to FIDO. 

In the Boyer-Lindquist coordinates V-/3 = 0, and under 
this condition equations (14,15,16,17) reduce to the corre- 
sponding equations in (Macdonald & Thorne 1982; Thorne 
et al.1986). For example, using equations (29,14) one finds 
that 

VxE = VxaE-£pB, (63) 
where 

Cf,& = {p-V)S-(&-V)p (64) 

is the Lie derivative of B along (3. Thus, equation (15) reads 

d t B - CpB + VxaE = 0, (65) 

which is equation (3.52) in Thorne et al.(1986). However, 
in other coordinate systems, e.g. the Kerr-Schild system, 
V-/3/0. 

The results presented in this section show that the 3+1 
equations of black hole electrodynamics can be written in 
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a more general and much simpler form than in (Thorne & 
Macdonald 1982; Macdonald & Thorne 1982). Our equations 
are remarkably similar to the familiar equations of electro- 
dynamics in matter. All effects of gravity are hidden in the 
constitutive equations 29,30) and in the spatial metric jij. 



3 STEADY-STATE FORCE-FREE 

MAGNETOSPHERES OF BLACK HOLES 

In this section we re-derive, within our 3+1 framework, some 
of the well known results for the steady-state force- free mag- 
netospheres of black holes obtained by Blandford & Zna- 
jek(1977). The only condition imposed on the coordinate 
system is d^g vll = dtg^v = and, thus, these results hold 
equally well both in the Boyer-Lindquist and the Kerr-Schild 
coordinates. 

The covariant force-free condition 

F„„J" = 

can be written in the 3+1 language as 

E J = 0, 

and 



(66) 



P E + JxB = 0, 
or 

D j = 0, 

and 

pD + j x B = 0. 

From eq.(69) it follows that 

D B = E B = 0, 



and 



3± 



DxB 



ExB 



(67) 



(68) 



(69) 



(70) 



(71) 



B 2 r B 2 ' 
where j ± is the component of electric current normal to B. 
Equation (71) states that this component is entirely due to 
the drift motion of charged particles. Since the drift velocity 
must be lower than the speed of light we have 

B 2 - D 2 > or B 2 -E 2 > 0. (72) 

Equations (70) and (72) are the 3+1 representations of the 
Lorentz invariant constraints of force-free approximation 

'F^F"" = and F^F» V > 

(Znajek 1977; Komissarov 2002a). Whereas the constraint 
(70) is automatically satisfied by any force-free steady-state 
solution, the constraint (72) is not and, hence, it always has 
to be checked (Znajek 1977). If this condition is not satisfied 
the Alfven wavespeed becomes complex and, thus, the sys- 
tem of force-free electrodynamics (magnetodynamics seems 
to be a better name) is no longer hyperbolic (Komissarov 
2002a). The physical reason why the force- free approxima- 
tion breaks down if E ■ B = but B 2 - E 2 < is quite 
obvious. Under such conditions, one can find a local iner- 
tial frame where the electromagnetic field is seen as a pure 
electric field. 

The conditions of axisymmetry and time-independency 
imply that 



Etf, — Ut,4> — U<t>,t = 0. 

Given this, equation (67) ensures that 

J p || Bp, E p _L Bp, 



(73) 



(74) 



where index p refers to the poloidal component of a vector. 
These show that there exists a purely azimuthal vector u> = 
Qdj, such that 



E — —lo x B 

and 



D = (w +(3)xB. 

a 



(75) 



(76) 



Substituting (75) into the stationary version of (15), one 
finds 



Vx(uxB) =0, 



which leads to 
B ■ Vfi = 0. 



(77) 



Thus, Q, which is called the angular velocity of magnetic 
field, is constant along the magnetic field lines. 

From equations (17) and (74) it immediately follows 

that 



B 6 d e Fl^ + B r d r H^ = 0, 



(78) 



and, thus, H$ is also constant along the magnetic field 
lines (This constant is denoted as Bt in Blandford & Zna- 
jek,1977.) 

Another parameter constant along the magnetic field 
lines is the scalar potential, which we denote as $. Indeed, 
from eqs. (39,75) one immediately finds that 



BV$ = 0. 



(79) 



Since the source terms in eqs. (50,51) vanish, the energy 
and the angular momentum of the electromagnetic field are 
conserved. 



dte + V-S = 0, 
dtl + V-L = 0. 



(80) 
(81) 



Substituting E from (75) into eq.(55) and eq.(54) one ob- 
tains the following expressions for the poloidal component 
of the angular momentum flux vector 



Lp HffyBp 



(82) 



and the poloidal component of the energy flux vector 

S p = -(H^Bp, (83) 

where B p is the poloidal component of the magnetic field. 
Thus, both the energy and the angular momentum are trans- 
ported along the poloidal field lines. 

Finally, from eqs. (29,75) one obtains the following use- 
ful result 



(B — D )a = B (a — (i ) +(w-B + f3-B) 2 

-B 2 {to 2 + 2w/3) . 



(84) 
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4 GENERALIZED OHM'S LAW 



It is generally accepted that the magnetospheres of black 
holes are filled with perfectly conducting e + -e~ plasma 
(Blandford & Znajek 1977; Phinney 1982; Beskin et al.1991; 
Hirotani & Okamoto 1998). Should, however, a current sheet 
be formed in the magnetosphere (see Komissarov, 2002) the 
perfect conductivity approximation would fail and a model 
of electric resistivity would be required. It is well known 
that magnetic field strongly effects the electric properties of 
plasma by suppressing electric conductivity across the mag- 
netic field lines and, thus, making the simple scalar Ohm's 
law inadequate. In fact, the total electric current splits into 
three components 



3 = °\\ E \\ + +3\ 



(85) 



where Eu is the component of electric field parallel to the 
magnetic field, E± is the component of electric field per- 
pendicular to the magnetic field, and j d is the drift current 
which is perpendicular both to the electric and to the mag- 
netic field vectors (Cowling 1976; Mestel 1999). The magne- 
tospheric plasma is collisionless and its resistivity arises via 
collective processes, "the anomalous resistivity", and emis- 
sion of photons, "the radiative resistivity". The radiative 
resistivity associated with the inverse Compton scattering 
of the photons emitted by the accretion disc has already 
been considered as a source of the internal resistivity in the 
potential gaps of the black hole magnetospheres (Beskin et 
al.1991; Hirotani & Okamoto 1998). Since, this is a rather 
simple and yet robust mechanism, it deserves to be consid- 
ered in details. 

If the gyration and the collisional time-scales are much 
smaller compared to the macroscopic time-scale then the 
mean velocity of a charged particle in electromagnetic and 
radiation fields, as observed by the local FIDO, is governed 
by 



q(E + vxB) - (W 2 



l)a T Ub — 
v 



0, 



(86) 



where the last term describes the inverse Compton interac- 
tion with soft background photons (Beskin et al.1991; Hi- 
rotani & Okamoto 1998). Here p, v, W, and q = ±e are the 
momentum, the velocity, the Lorentz factor, and the elec- 
tric charge of the particle respectively, ot is the Thomson 
cross section, and u b is the energy density of the radiation 
field. For simplicity, we assume that the radiation field is 
isotropic (Since the Boyer-Lindquist FIDO becomes singu- 
lar near the event horizon this condition cannot be satisfied 
there at any length. However, for the Kerr-Schild FIDO this 
seems to be more or less acceptable everywhere.) One can 
solve this equation for the perpendicular and the parallel to 
the magnetic field components of velocity to obtain 



±-E n 

X 

, e 



■E± + 



ExB 



1 + S 2 B 2 



where 



X = mWvc, 
8 = v b /v c , 
where 



(87) 
(88) 

(89) 
(90) 
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Figure 1. Schwarzschild black hole in uniform magnetic field. 
The black circle shows the position of the event horizon and the 
continuous lines show the magnetic flux surfaces. 



eB 



(91) 



(92) 



mWc 

is the electron gyration frequency, 

_ (W 2 - l)a T u b 
mWv 

is the effective frequency of the inverse Compton collisions; 
here m is the particle rest mass. The last term in eq.(88) does 
not depend on the sign of the electric charge and describes 
the usual drift motion effected by the collisions. 

Equations (87,88) give us the generalized Ohm law of 
the same form as eq.(85) with 



ne 



o± = 



3d = P 
and 



mWvc 
1 

TTs 2 ' 

8 2 



ExB 



1 + S 2 B 2 



p = e(n + — n_), n = n + +n_, 



(93) 
(94) 
(95) 

(96) 



where n± are the number densities of electrons and 
positrons. 

The typical values of B and u b in vicinity of a su- 
permassive black of M w 10 8 M Q are B w 10 4 Gauss and 
u b ~ 10 7 ergs _1 . Hence, we have 

-i 



0.1W,s" 



Vb ~ 10 W S , 



and 
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10 12 W~ 



(97) 



where W is the typical Lorentz factor of charged particles. 
Thus, the conductivity across the magnetic field is strongly 
suppressed indeed, and the drift velocity is given by the 
familiar expression 

v d = ExB/B 2 . 

For the Goldreich- Julian number density, uqj w 0.1 cm~ 3 , 
we obtain 

<T|| = W a W' 2 s- 1 a ± = HT 16 W"V\ 
whereas the macroscopic frequency is 



lrrV 1 . 



2MG 



Thus, one may confidently put E« = and ignore the cross- 
field conductivity current, a±E± = 0. This corresponds to 
the limit of force-free electrodynamics. Indeed, under these 
conditions 

pE + jxB = P E ± + j d x B = 0. 

In this limit, the role of the E± component of electric field 
is reduced to driving of the drift current. 

Inside a current sheet the cross-field conductivity can no 
longer be ignored. In fact, it has got to be governed by a self- 
regulatory mechanism ensuring marginal screening of the 
electric field. Indeed, let us assume, for a moment, that the 
unscreened component of electric field is of the same order 
as the magnetic field. Then from equation (86) it follows 
that 



W i 



eB 
arUb 



10° 



(98) 



This implies 8 « 1 (see eq.97 ) and, thus, similar conductivi- 
ties both along and across the magnetic field. However, par- 
ticles with such a high Lorentz factor are bound to launch co- 
pious pair cascades (Beskin et al.1991; Hirotani & Okamoto 
1998). This would dramatically increase the electric con- 
ductivity and lead to effective reduction of the unscreened 
component of the electric field. Thus, the electromagnetic 
field of a magnetospheric current sheet is expected to be 
very close to the state of marginal screening with 

B 2 — E 2 = 0. (99) 

(In fact, the electric field should remain slightly stronger 
than the magnetic field.) Since in this state, the typical par- 
ticle's Lorentz factor is much smaller than the one given by 
cq.(98), one should have 



cm > fx, E± 3> -Ell 



and 



3d = P- 



ExB 



B 2 



(100) 



Summarizing, the cross-field conductivity is essentially 
zero in the main body of the black hole magnetosphere, 
where it is basically force-free, and dramatically increases 
in current sheets to ensure that the magnitude of £i ex- 
ceeds the magnitude of the magnetic field only by a small 
margin. 
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Figure 2. Monopolc field solution for a = 0.1 and Bo = 1 at time 
t = 50. Panel a) shows the angular velocity of magnetic field lines 
at r=3. The perturbative solution of Blandford and Znajck gives 
f! = 0.5f2f,, where Clh is the angular velocity of the black hole. 
Panel b) shows of the numerical solution (the crosses) and of 
the perturbative solution (the continuous line) at r = 3. Panel c) 
shows E ■ B along 8 = 1 and 9 = n — 1 . The small unscreened 
component of electric field drives the conductivity current 
towards the black hole in the upper hemisphere and away from it 
in the lower hemisphere. 
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5 NUMERICAL STUDY OF BLACK HOLE 
MAGNETOSPHERES 

To describe the black hole space-time we adopted the Kerr- 
Schild coordinates, {t, <f>,r,6}, which are free from the coor- 
dinate singularity at the event horizon hampering the Boyer- 
Lindquist coordinates (see Appendix B). This allows us to 
put the inner boundary of the computational domain inside 
the horizon and impose radiative boundary conditions on 
this boundary. 

All simulations described below are axisymmetric. The 
computational grid is uniform in the ^-direction with A9 = 
n/ne, where ng is the total number of cells in this directions. 
The cell size in the r-direction, Ar, is determined via the 
condition of the same length in both directions, 7 rr Ar 2 = 
7ee A6 2 . The time step is determined by the Currant stability 
condition. 

The details of the resistivity model used in these simu- 
lations are given in Appendix C3. 

5.1 Wald's solution 

Wald (1974) obtained the following vacuum solution for a 
rotating black hole immersed into a uniform magnetic field 
aligned with hole's rotational axis: 

= B (m M + 2ak M ), (101) 

where k" — d t and rrf = 84, are the Killing vectors of the 
Kerr spacetime. 

In the case of a Schwarzschild black hole eq.(lOl) reads 

E = 0, B = ^p(0,-7 W , fl ,7 W , r ). (102) 

Outside of the horizon this solution satisfies the plasma 
equations (14-17,85) with zero p and j. Indeed, from eq.(29) 
one finds 

D = --(3xB. 

a 

Thus, Du = and no electric current is driven along the 
magnetic field. Provided B 2 — D 2 > no electric current 
will be driven across the magnetic field lines either. From 
eq.(84) one has 




which is strictly positive outside of the event horizon. (Inside 
of the horizon a 2 - 2 < and B 2 - D 2 becomes negative 
near the equatorial plane where B r = 0. However, these 
details do not effect the exterior solution.) 

The computational grid in this test problem has 100 
cells in the ^-direction and 80 cells in the r-direction covering 
the domain r € [1.4, 29]. The initial solution is described by 
cq.(102). Figure 1 shows the magnetic flux surfaces of the 
numerical solution at time t = 5 as well as the flux surfaces of 
the exact steady-state solution (102). In fact, these solutions 
are so close that one cannot see the difference between them 
in this figure. 

5.2 Blandford-Znajek's monopole solution. 

Using a perturbation method Blandford and Znajek(1977) 
constructed a perturbative analytic solution for a slowly ro- 



tating black hole with monopole magnetic field that matches 
the flat space solution of Michel(1973) at infinity and sat- 
isfies Znajek's "boundary condition" (see Sec. 6.2) imposed 
at the event horizon. The zero order solution is the one for 
a nonrotating black hole with purely radial magnetic field 

B r = B sin 0/^/7. (103) 

If a < 1 then the poloidal magnetic field is not expected 
to be very different from purely radial with B r given by 
(103) and, thus, B r ^fy remains constant along the field line. 
Moreover, the Znajek condition (eq.U8) reads 

= (n-n h )sm6B r ^j, (104) 

where Q,h = a/(r+ + a, 2 ) w a/4 is the angular velocity of the 
black hole. On the other hand, Michel's monopole solution 
in flat spacetime gives 

H 4> = -VLsm6B r ^. (105) 

Matching of given by (104) with H$ given by (104) 
leads to 

n = 0.5fi h = a/8 and = sin 2 9 (106) 

8 

In these simulations case we used Bq = 1 and put 
a = 0.1 because for this value of a the corresponding force- 
free numerical solution is known to be very close to the per- 
turbative one (Komissarov 2001b). The computational grid 
in this test problem had 100 cells in the ^-direction and 150 
cells in the r-direction covering the domain r £ [1.4,260]. 
We tried two very different initial solutions for this problem. 
One of them describes the vacuum monopole field (this solu- 
tion was found numerically.) The other one has the same B p 
as in the Blandford-Znajek solution and E = 0. However, 
in both these cases the outcome was the same, the numeri- 
cal solution gradually evolved towards the Blandford-Znajek 
one. Figure 2 shows the angular distribution of £1 and H$ 
at r — 3 as well as the distribution of the parallel compo- 
nent of the electric field in the northern and in the southern 
hemi-spheres at t — 50. 

5.3 Split-monopole solution 

Since the magnetic charges do not seem to exist, the 
monopole solution of Blandford and Znajek (1977) is rather 
artificial. A small modification, namely alternation of the 
magnetic field direction in one of the hemi-spheres, seems 
to overcome this problem. The result, known as the "split- 
monopole" solution, is still a force-free solution but now with 
zero total magnetic flux. However, such a configuration is 
meaningful only if there is a massive perfectly conducting 
thin disc in the equatorial plane of the black hole. Other- 
wise, the equatorial current sheet is not sustainable - the 
magnetic field is expected to reconnect and "fly away" . In- 
deed, this is what is observed in numerical simulations (see 
Figure 3.) In this problem a = 0.1 and we use almost the 
same grid as in the Wald problem (Sec. 5.1). 

5.4 The Magnetospheric Wald Problem for 
rotating black holes 

The Blandford-Znajek solution is the only global analytical 
solution for magnetospheres of rotating black holes found 
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Figure 3. Escape of the split-monopolc magnetic field from a Schwarzschild black hole. Left panel: Magnetic flux surfaces of the split- 
monopole solution which was used as an initial solution in these numerical simulations. Right panel: Magnetic flux surfaces of the 
numerical solution at t = 5. 



so far and for this reason it has been playing a key role in 
the development of the black hole electrodynamics. One im- 
portant property of this solution is that all magnetic field 
lines penetrate the black hole horizon. Macdonald (1984) at- 
tempted to construct numerical steady-state solutions for a 
more reasonable configuration of magnetic field where only 
a fraction of magnetic field lines originates from the black 
hole itself. The remaining magnetic flux splits between the 
field lines originating from the accretion disk and the field 
lines passing through the gap between the hole and the disc. 
In general, the angular velocity of magnetic field lines in 
steady-state force-free magnetospheres has to be prescribed, 
so one faces the task of setting physically sensible boundary 
conditions for all these three different types of magnetic field 
lines. In the case of the field lines originating from the ac- 
cretion disc the solution is obvious. Their angular velocity is 
given by the angular velocity of the disc at the foot points. 
As for the other two kinds of magnetic field lines, this task is 
less trivial. In their solution, Macdonald and Thorne(1982) 
and later Macdonald (1984) appealed to the existing analogy 
between the black hole horizon and a rotating conducting 
sphere. They concluded that only the field lines penetrating 
the event horizon rotate, whereas in the gap f2 = 0. 

A somewhat simpler problem is the magnetospheric 
(plasma-filled) version of the Wald (1974) problem for a ro- 
tating black hole (see also Sec. 5.1). In this problem, just 
like in the problem considered by Macdonald (1984), only 
a small fraction of magnetic field lines penetrate the black 
hole horizon. If the analysis of Macdonald & Thorne (1982) 



was correct then only these field lines would be forced to ro- 
tate. Komissarov(2002b) tried to find a steady state force- 
free solution to this problem by means of time- dependent 
numerical simulations but failed. The numerical solution in- 
variably evolved towards the state where B 2 — D 2 turned 
negative inside the ergosphere. In fact, the solution seemed 
to indicate the development of a current sheet in the equato- 
rial plane within the ergosphere with all magnetic field lines 
penetrating the current sheet being forced to rotate in the 
same sense as the black hole. If this conclusion is correct, a 
critical revision of the current perception of the role of the 
event horizon in the black hole electrodynamics, as well as of 
the virtues of the Membrane Paradigm, is required. Thus, 
the magnetospheric Wald problem is an ultimate Rosetta 
Stone for the research into the black hole electrodynamics. 

To achieve high resolution within the ergosphere these 
simulations were carried out using the multi-grid technique. 
We start with a relatively low resolution grid and continue 
simulations until the solution becomes more or less steady 
within r = 4. Then the resolution is increased by the factor 
of 2 and the simulations are continued until new approxi- 
mately steady-state solution is reached and so on. During 
the grid refinement the numerical solution on the finer grid 
is found via interpolation. The final grid has 800 cells in 
the ^-direction (9 G [0, 7r]) and 1000 cells in the r-direction 
(r G [0.9r+, 110]). The initial solution is described by the 
same B as in the vacuum solution of Wald (1974, eq.101) 
and has E — 0, which implies a non-rotating magnetosphere. 

Figure 4 shows the final solution, at t = 126, for a 
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Figure 4. Magnetospheric Wald problem. Left panel: The angular velocity of magnetic field lines. There are 15 contours equally spaced 
between and 0.67. The angular velocity first gradually increased towards the axis but then reaches a maximum and goes slightly down. 
The thick lines show the ergosphere (the outer line) and the inner light surface (the inner line). Middle panel: The magnetic flux surfaces. 
Right panel: The distribution of (B 2 - D 2 )/max(B 2 , D 2 ). There are 15 contour equally spaced between -0.12 and 1.0. This quantity 
monotonically decreases towards the current sheet in the equatorial plane within the ergosphere. The thick line shows the ergosphere. 



Kerr black hole with a = 0.9. As suggested in Komis- 
sarov(2002b), a current sheet is formed in the equatorial 
plane within the black hole ergosphere. This is clearly seen 
in the right panel of fig. 4 which shows the distribution of 
(B 2 — D 2 )/max(B 2 , D 2 ). Near the equator the predomi- 
nantly radial electric field is larger than the magnetic field 
and drives the electric current across the poloidal magnetic 
field lines. Both the radial component (the middle panel of 
fig. 4) of magnetic field and its azimuthal component exhibit 
a break in the equatorial plane on the scale of the current 
sheet. The most important result is shown in the left panel 
of fig. 4: all magnetic field lines penetrating the ergosphere 
are forced into rotation in the same sense as the black hole 
irrespective of whether they eventually cross the event hori- 
zon or not. Along these field lines there are outgoing fluxes 
of both energy and angular momentum. Indeed, in steady 
state force-free magnetospheres the angular momentum flux 
is proportional to and the energy flux is proportional to 
QH,/, ( see eqs.82,83. Notice that Q and H^, are the same in 
the Boyer-Lindquist as in the Kerr-Schild coordinates; see 
Appendix B). As one can see in fig. 5, both these quantities 
are nonvanishing along the field lines penetrating the ergo- 
sphere. Within the current sheet the electromagnetic energy 
and angular momentum are not conserved and, thus, the 
numerical results suggest that it is the current sheet that 
supplies both the energy and the angular momentum for 
the force-free magnetosphere above and below the sheet. 



To show that this makes sense let us consider the 
sources of energy and angular momentum in the symmetry 
plane, 8 — tv/2. Because of the symmetry D\\ vanishes in this 
plane. Thus, we can ignore the dissipation due to au , and 
we can still introduce vector uj to describe D as in eq.(76). 
In addition, the marginal screening of electric field implies 
B 2 — D 2 . The source terms in the Boyer-Lindquist coordi- 
nates have the same sign as in the Kerr-Schild coordinates 
and this allows us to utilize the Boyer-Lindquist coordinates 
to simplify calculations. 

The source term in the energy equation (50) is 

m = -e j. 

Substituting the expressions for E and J from (29,62,85) 
into this equation one obtains 



-a 2 a ± D 2 - aa ± D-((3 x B) 



Then, provided D. < D. z , where Sl z = —g^t /g<t>4, is the angular 
velocity of the "zero angular momentum observer" (ZAMO, 
Bardeen et al.,1973), we have 



D-((3xB) = -\(3\B 2 . 



Hence, we obtain 



W t = aa ± B 2 (\l3\-a). 



(107) 
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Figure 5. The distribution of the angular velocity of magnetic lines (left panel) and (right panel) in the magnetospheric Wald 
problem at r = 1.8. 



Similar calculations lead to the expression for the angular 
momentum source, W^,, in the equatorial plane 

W4, = -(JxB)-m = aa ± B 2 . (108) 

Thus, provided the magnetosphere rotates slower then 
ZAMO, the resistive dissipation in the ergospheric current 
sheet can provide positive sources of electromagnetic energy 
and angular momentum. Since, the inner light surface en- 
closes the current sheet (see the left panel of figure 4) the 
condition Q < Q z is satisfied in our numerical solutions. Ob- 
viously, this source of electromagnetic energy and momen- 
tum implies the existence of the same magnitude sink of en- 
ergy and angular momentum for the radiation field. Thus, 
the energy transfer has to be essentially the same as pro- 
posed in (Punsly 2001). Namely, the electromagnetic field 
pushes plasma particles into orbits with negative energy at 
infinity, these particles push background photons into orbits 
with negative energy at infinity and finally these photons 
are swallowed by the black hole. 

As the grid resolution increases the current sheet be- 
gins to show sings of unsteady behavior. Whether this is 
only a numerical effect or an indication of the current sheets 
instability remains to be investigated. 



6 DISCUSSION. THE NATURE OF 

BLANDFORD-ZNAJEK MECHANISM 

The results of numerical simulations presented in Sec. 5 allow 
us to draw the following two general conclusions concerning 
the electrodynamics of black hole magnetospheres. 

First of all, the rotational energy of black holes can 
indeed be extracted electromagnetically. In particular, the 
Blandford and Znajek solution is certainly stable and, hence, 
does not clash with causality. 

On the other hand, there is much more to the electro- 
dynamics of black holes than it is proposed in the Mem- 
brane Paradigm and the analogy between a magnetized ro- 
tating conducting sphere and the black hole horizon is at 



least incomplete. In fact, closer inspection of the causality 
arguments due to Punsly and Coroniti shows that it is the 
Membrane Paradigm which is most directly under attack. 
The Blandford-Znajek solution itself is involved mainly be- 
cause it is widely considered as inseparably linked with the 
paradigm. 

6.1 The origin of torque 

6.1.1 Gravitationally induced electric field. 

It is well known that a black hole with zero total electric 
and magnetic charge cannot have its own magnetic field 
and, ultimately, any magnetic field penetrating the hole's 
ergosphere has to be supported by external currents. This 
fact alone makes black holes very different from magnetized 
stars like pulsars. Moreover, a steady state axisymmetric 
vacuum electromagnetic field cannot be used to extract en- 
ergy and angular momentum of a rotating black hole. In- 
deed, the steady-state vacuum equations, 

VxH = 0, VxB = 0, 

ensure 

H 4 , = E 4 , = (109) 

for axisymmetric configurations. Then from eqs. (54,55) one 
finds that the poloidal components of the energy and the 
angular momentum flux vectors vanish as well. In order to 
extract energy and angular momentum the electromagnetic 
field has to be modified by magnetospheric charges and cur- 
rents. However, in order to drive such currents and, perhaps, 
even to create charged particles via pair cascade in the first 
place (Blandford & Znajek 1977; Beskin et al.1991), the elec- 
tric field should not be screened. That is at least one of the 
conditions (70,72) has to broken in the vacuum solution. As 
it is well known, the vacuum solution due to Wald (1974) 
has this property but we need to know whether it is generic. 

Theorem. There are no steady-state axisymmetric vac- 
uum solutions supported by remote sources of magnetic field 
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that simultaneously satisfy both 

D B = and B 2 - D 2 > 

along the magnetic field lines penetrating the ergosphere of 
a rotating black hole. 

To prove this statement let us assume that electric cur- 
rents of these sources can be arranged in such a way that 
DB = everywhere in the vacuum zone surrounding the 
black hole. Just like in the case of a force- free magneto- 
sphere this condition ensures the existence of the vector field 
w = £lm such that O is constant the along magnetic field 
lines and 

E = -ljxB. 

The remote sources may be regarded as non-rotating relative 
to the black hole and, thus, E = along the field lines 
penetrating the ergosphere. Then from eq.(29) one obtains 

D= -(3xB. (110) 

a 

Thus, even though there are no external sources of electric 
field, the electric field is generated near the black hole as 
the result of the inertial frame dragging effect manifested 
in the shift vector [3. This is not entirely unexpected. Rela- 
tive to any physical observer located within the ergosphere, 
including FIDO, these external sources do rotate. In the 
Boyer-Lindquist coordinates f3 is a purely azimuthal vector 
and B,j, = H^/a — (see eq.109). This allows to find 

B 2 -D 2 = ?—^B 2 . (Ill) 

Ol 

Within the ergosphere a 2 — (i 2 < and, hence, B 2 — D 2 < 0. 
Thus, no matter how strong is the magnetic field within 
the ergosphere, the gravitationally induced electric field is 
always stronger. This completes the proof. 



Obviously, these external sources may not be infinitely 
far away and strictly non-rotating. In such a case, Q is not 
exactly zero and instead of eq.(lll) one obtains 

a 2 (B 2 -D 2 ) = ~B 2 f(Q,r,e), (112) 

where 

/(fi, r, 6) = g^Q 2 + 2g t4 ,Q, + gtt 

is the light surface function (see Appendix A) . Thus, the con- 
dition B 2 — D 2 > is equivalent to the condition Q > Q m in 
(see eq.A9). Since within the ergosphere Q. m in > and 
reaches the value of Q h at the event horizon, this requires 
the remote sources to rotate really fast. In fact, unless the 
remote sources of magnetic field are located within the outer 
light surface corresponding to f2 m i«, which may be theoreti- 
cally possible but which is highly unlikely to occur naturally, 
this condition cannot be satisfied. 

Thus, the vacuum field created within a black ergo- 
sphere by distant sources must have unscreened electric field 
capable of driving electric currents, provided the charged 
particles are injected there somehow. The most popular 
mechanism involves an e _ -e + pair cascade in strong elec- 
tric and radiation field (Blandford & Znajek 1977; Beskin et 
al.1991). Let us consider an initial vacuum solution of the 
kind discussed above, that is a solution with D B = and 
= 0, and figure out what occurs when plasma is injected 
(In fact, our initial numerical solutions in the Bladford- 
Znajek problem and the magnetospheric Wald problem are 
of this type.) 

First of all, FIDO's electric field E = D given by (110) 
will drive the conductivity current across the magnetic field 
lines within the ergosphere. This will result in the electric 
charge separation and the drop of the electrostatic potential 
along the magnetic field lines entering the ergosphere (see 
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figure 6). The induced parallel component of the electric 
field will drive the conductivity current along the magnetic 
field lines leading to the dynamo of i/^ according to the 
Ampere law V x H — J. As one can see in the left panel 
of figure 6, the sign of the generated is exactly the one 
which is required to slow down the black hole. However, 
we need to verify that this is not a temporary phenomenon 
and a nonvanishing and, hence, a nonvanishing poloidal 
current, will be a property of the magnetosphere when it 
relaxes to a steady-state. 

6.1.2 Poloidal currents in static magneto spheres. 

Let us suppose that the magnetic field lines penetrating the 
ergosphere are anchored to distant plasma of infinitely large 
inertia. Then the created magnetosphere will eventually be- 
come static with SI = 0, E = 0. From eq.(84) one finds that 
in a force-free region of such magnetosphere 

B 2 -D 2 =B l C? ~f +\(f3.Bf. (113) 

a 2 ex 

In the Boyer-Lindquist coordinates this reads 
a 2 (B 2 -D 2 ) = B 2 (a 2 -/3 2 )+(^pj Hi (114) 

The difference between this result and equation (111) is en- 
tirely due to the azimuthal magnetic field created by the 
poloidal electric current. According to equation (76), the 
growth of this component of magnetic field does not lead to 
the growth of D. 

A number of important conclusions can be drawn from 
this simple result. 

• First of all, inside the ergosphere a 2 — ft 2 < and for 
B 2 — D 2 to be positive must be nonvanishing. This shows 
that screening of electric field by the injected plasma within 
the ergosphere implies generation of poloidal electric cur- 
rents. These currents give rise to a nonvanishing and, 
thus, to a nonvanishing flux of angular momentum. Zero 
angular velocity of the magnetosphere means zero Poynting 
flux from the hole (see eq.83). Thus, the total mass of the 
hole, M, remains constant but its irreducible mass, Mi rr , 
grows due to the loss of angular momentum. One may com- 
pare this case with a shorted out Faraday disk where its 
rotational energy is being used only for Ohmic heating of 
the disc itself. 

• Secondly, it makes no difference whether the magnetic 
field lines penetrate the event horizon or not, the only fac- 
tor that matters is whether they penetrate the ergosphere. 
The analogy between the event horizon and the rotating 
conducting sphere emphasized in the Membrane Paradigm 
does not reveal this important fact and, therefore, proves to 
be of rather limited value. 

• Finally, it explains the origin of the equatorial current 
sheet in the magnetospheric Wald problem. Indeed, the mag- 
netosphere is symmetric relative to the equatorial plane, and 
for this reason alone must vanish in this plane. Thus, 
the B 2 — D 2 > condition cannot be satisfied in the ergo- 
spheric part of the plane and the force-free approximation 
breaks down. As we have seen in Sec. 5. 4 , this current sheet 
serves as a localized source of angular momentum for the 
surrounding force-free magnetosphere. Moreover, it provides 



the required closure of the electric circuit at the black hole 
end (see the right panel of fig. 6). 

6.1.3 Poloidal currents in rotating magneto spheres 

Since the assumption of infinite inertia of the surrounding 
plasma, made in Sec. 6. 1.2, is unrealistic and was made only 
for the sake of argument, the outgoing flux of angular mo- 
mentum will inevitably result in rotation and, thus, extrac- 
tion of the black hole energy as well. The equilibrium value 
of Q depends of the details of the interaction. In our simu- 
lations, the black hole was surrounded by massless plasma 
and the equilibrium value of fi was determined by the rate 
of deposition of energy and angular momentum in the sur- 
rounding space by means of propagating waves. 

Even after reaching a force-free equilibrium the magne- 
tospheric electric field does not become completely screened. 
Within the ergospheric current sheet, E± always remains 
slightly stronger than the magnetic field and keeps driving 
the cross field conductivity currents, thus, sustaining the po- 
tential drop along the magnetic field lines entering the cur- 
rent sheet. In the force- free region itself it is the small resid- 
ual component of E\\ that drives the poloidal conductivity 
currents (see fig.2c), just as it occurs in the wires connected 
to the Faraday disc. 

It is also quite clear why the magnetic field lines remain- 
ing outside of the ergosphere along their entire length do 
not rotate. Although the vacuum Wald solution has an un- 
screened E\\ component of the electric field which can trigger 
pair production and drive poloidal currents for a while, the 
injected charges eventually redistribute along the field lines 
and screen E\\. Then, because E± is too weak (see eq.114), 
even in the equatorial plane, to drive the cross field cur- 
rent, the poloidal currents die out completely. The force free 
conditions are met even when H$ — 0. 

The spatial grids of both the Boyer-Lindquist and the 
Kerr-Schild coordinates, though stationary relative to dis- 
tant inertial observers at rest relative to the black hole, ro- 
tate superluminally relative to any local physical observer 
located within the ergosphere. This property of black hole 
space-time is crucial in determining the properties of their 
magnetospheres. To illustrate the point let us consider a 
steady-state magnetosphere rotating with angular velocity 
f2 in flat space-time (we assume that some particle injection 
mechanism ensures plentiful supply of electrical charges re- 
quired to screen the electric field.) Let us introduce a spher- 
ical grid rotating with the same angular velocity. The cor- 
responding coordinate transformation 

t = t', r = r',4> = <p'- fit', 6 = 9', 

where the primed coordinates refer to the nonrotating grid, 
results in the metric form 

ds 2 = (/3 2 -a 2 )dt 2 + 2f3 4> d<f>dt + dr 2 + r 2 (sin 2 8 d(f> 2 + d(9 2 )(115) 
where 

a = 1, (3 2 — f2 2 sm 2 #r 2 , = Q.sin 2 6r 2 . 

We notice that this metric is rather similar to the Kerr met- 
ric as the surface rsin# = vj c — 1/Q separates the inner 
region where the time coordinate is time-like from the outer 
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region where it is space-like (obviously, uj c is the cylindrical 
radius of the light cylinder.) Relative to the rotating grid 
the magnetosphere is static and, thus, E = 0. Calculations 
similar to those leading to eq.(114) give us 

B 2 -D 2 =B 2 {l-{vj/vj c ) 2 ) + Q. 2 Hl. (116) 

Thus, high conductivity ensures — 7^ outside of 
the light cylinder. The similarity with black hole magneto- 
spheres is startling. 

The result (116) serves another goal. It shows that H<f, 
must remain nonvanishing in rotating magnetospheres of 
black holes for any value of / 0. 

6.I.4 Where is the unipolar inductor of a Kerr black 
hole? 

In the case of the Faraday disc, which is a classical example 
of a unipolar inductor, electrons are forced to participate 
in the disc rotation via collisions with other disc particles. 
This provides the electromotive force, qvxB, which results 
in electric charge separation and, hence, the voltage drop 
between the disc rim and its centre. When the disc is used 
as battery in a closed electric circuit the electromotive force 
continues to push electrons against the electric force and 
across the magnetic field lines. This is essential for sustaining 
the potential drop and for providing the current closure. 
Spinning magnetized cosmic objects like stars or accretion 
discs generate an electric field in the very much same way. 
Such a field is often described as rotationally induced. 

Although, the Membrane paradigm invites to treat 
black holes in a similar fashion, in reality they are rather 
different. Both in the Blandford-Znajek solution and in our 
solution to the magnetospheric Wald problem there is no 
any massive conducting rotating object and, thus, no usual 
electromotive force driving electric current over its surface. 
In spite of this curious feature the potential drop across the 
magnetic field lines still exists and the electric currents still 
flow. 

The answer to this paradox seems to reside in another 
peculiar property of rotating black holes - contrary to our 
everyday experience the electric charge separation in not the 
only way of creating stationary electric field in their vicinity. 
The vacuum solution found by Wald(1974) shows that such 
an electric field can be induced gravitationally. Thus, the 
usual electromotive force is not required at least for initial 
charging of the black hole "battery" . Punsly and Coronity 
(1990a) argued that such electric field can not sustain sta- 
tionary electric currents. Their first point is based on the 
change of sign of the parallel component of the electric field 
in the exact Wald solution. Indeed, this seems to be incon- 
sistent with the electric current flowing along the poloidal 
magnetic field line in one particular direction. However, the 
electric field modified by the magnetospheric charges does 
not have to retain this property. In their second point, they 
argue that the mere fact of existence of any initial electric 
field does not matter very much. For example, a capacitor 
can drive only a transient electric current. Unless a battery is 
continuously recharged via the action of some electromotive 
force it cannot drive a stationary electric current. In other 
words, the initial transient currents can redistribute elec- 
tric charge in the black hole magnetosphere in such a way 



that the electric field becomes completely screened and can 
no longer drive the electric currents. However, as we have 
shown above, such a final state is not possible within the 
black hole ergosphere. Black holes do not allow stationary 
solutions with screened electric field and vanishing poloidal 
electric currents. 

Thus, the electrodynamics of rotating black holes is 
very different from electrodynamics of usual magnetic ro- 
tators and their batteries operate on other principles than 
the classical unipolar inductor. The key role played by the 
ergosphere in black hole electrodynamics allows us to call it 
the "driving force" of the Blandford-Znajek mechanism. 

6.2 Causality arguments 

Apart from the absence of the proper unipolar inductor Pun- 
sly and Coronity (Punsly & Coroniti 1990a; Punsly 2001) 
criticized the electromagnetic mechanism on the ground of 
causality principle. The key point of their criticism is con- 
cerned with the role of Znajek's horizon condition in the 
Blandford-Znajek solution (Blandford & Znajek 1977, Zna- 
jek 1977; see also Sec. 5. 2.) Since Blandford and Znajek em- 
ployed the Boyer-Lindquist coordinate system, the horizon 
appears as a singular boundary of their spatial domain and 
the Znajek condition appears as a boundary condition. In 
fact, this perception of the Znajek's condition is fully acco- 
modated in the Membrane paradigm, where this condition is 
used to endow the event horizon with the properties of rotat- 
ing conducting surfaces, reintroducing the missing unipolar 
inductor in a somewhat ghostly form. 

Punsly and Coroniti (1990a) convincingly argued that 
this analogy must be wrong. They pointed out that the BZ- 
solution describes not just an outgoing wind, as in the case 
of pulsar winds, but also an ingoing wind that passes though 
its own critical surfaces (see also Takahashi et al. 1990 and 
Appendix A.) As the result, the event horizon is causally 
disconnected from the outgoing wind and, thus, not suit- 
able for imposing boundary conditions that would determine 
the properties of the outgoing wind. This made Punsly and 
Coroniti (1990a) to conclude that although the Blandford- 
Znajek solution may be a valid solution of steady-state equa- 
tions it must be unstable and thus the time-dependent solu- 
tion would never evolve towards it. However, the results of 
numerical simulation presented in (Komissarov 2001b) and 
in this paper show just the opposite - the Blandford-Znajek 
solution is, in fact, asymptotically stable and, hence, causal. 

The answer to this controversy seems to be very simple. 
Znajek's horizon condition is not a boundary condition af- 
ter all. Indeed, Znajek derived the horizon condition from a 
very natural requirement. Namely, the electromagnetic field 
at the event horizon had to be non-singular when measured 
by a local free-falling observer. However, in the limit of force- 
free electrodynamics, the event horizon coincides with the 
fast critical surface of the ingoing wind. This immediately 
follows from the fact that the fast wave of force-free electro- 
dynamics propagates with the speed of light (Komissarov 
2002a). All this suggests that Znajek's condition is in fact 
a regularity condition at the fast critical point of the ingo- 
ing wind, obtained by Znajek in a somewhat unconventional 
way. 

There is an important difference between a boundary 
condition and a regularity condition. Boundary conditions 
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are set on the boundaries of computational domain and 
select particular solutions to both steady-state and time- 
dependent problems. Regularity conditions apply only to 
steady-state problems when one is looking for solutions pass- 
ing smoothly through critical points, at which the steady- 
state equations change their type. 

In the Kerr-Schild coordinates, where there in no co- 
ordinate singulatity at the horizon, the critical nature of 
the Znajck condition becomes very clear. Let us consider 
a steady state force-free solution smoothly passing through 
the event horizon. In such a solution the angular velocity of 
magnetic field, f2, and are constant along the magnetic 
field lines. Following Weber & Davis(1967) we can use these 
constants in order to find the relationship between and 
B r along a given magnetic field line. First, we eliminate D 
from eq.(30) using eq.(76) to obtain 

H = aB + 1 [uifi-B) - B(J3-u + 2 ) + (3(f3-B)] . 

From this equation we derive, after rather involved calcula- 
tions, the following result: 



B 



aH^ - B r sin 2 9(2rtt - a) 



(117) 



A sin 2 9 

The denominator of the right hand side of this equation van- 
ishes at the event horizon, where A = 0. For B^ to remain 
finite the numerator has to vanish as well and we obtain 

Tr sin 2 9 _ (2r + fl — a) sin# 

Hs = (2r+n~a)B r = 2 , , ' a Am, (U8) 



a+ 



+ a z cos 2 



which is exactly the Znajek condition (Notice that H$, A$, 
and f2 are the same in the Kerr-Schild coordinates as in the 
Boyer-Lindquist coordinates; see Appendix B) Thus, Zna- 
jek's "boundary condition" is indeed a regularity condition 
at the fast critical point of the ingoing wind. In MHD ap- 
proximation, in the case where the particle inertia is not 
negligibly small, the fast surface is located outside of the 
event horizon and Znajek's horizon condition is no longer 
needed to determine the steady-state wind solution (Beskin 
& Kusnctsova 2000). 

In time, the dispute over causality of the Blandford- 
Znajck solution, and causality of the electromagnetic mech- 
anism as such, gradually evolved into a dispute over the 
role played by different waves in the development of the 
system of poloidal currents in black hole magnetospheres, 
e.g. (Blandford 1991; Punsly 1996; Punsly 2001; Blandford 
2002). The origin of this slight change of subject resides in 
the fact that in the Membrane paradigm the inner boundary 
conditions are imposed not at the actual event horizon but 
at the so-called "stretched horizon" , located just outside of 
the actual one (Thorne et al.1986). This allows communica- 
tion between the stretched horizon and the outgoing wind 
by means of fast waves giving hope of salvaging the per- 
ception of the event horizon as a unipolar inductor of black 
holes (Blandford 1991; Blandford 2002). Punsly (1996; 2001) 
dismissed this idea showing that in order to provide the nec- 
essary adjustment of poloidal currents the stretched event 
horizon would have to communicate with the outgoing wind 
by means of Alfven waves. However, this is not possible as 
the inner critical Alfven surface in the Blandford-Znajek so- 
lution is located well outside of the stretched horizon. 

The conclusion of Punsly (2001) on the crucial role of 
Alfven waves in establishing the global current system seems 



to be correct in general. This is particularly easy to see in the 
limit of force-free electrodynamics. In this approximation 
the electric charge and the electric current density can only 
be changed by hyperbolic waves. The fast waves cannot do 
this as they have identical properties to linearly polarised 
waves of vacuum electrodynamics (Komissarov 2002a). The 
only other option is the Alfven waves. 

However, once we have seen that 1) the event horizon 
does not play the role of a unipolar inductor, 2) the Znajek's 
horizon condition is just the usual regularity condition, and 
3) the key role in the electrodynamic mechanism is played by 
the black hole ergosphere, this dispute has to be somewhat 
redirected. What we need to verify is that the ergosphere is 
causally connected with the outgoing wind. In fact, as it is 
shown in Appendix Al, the inner Alfven surface is always 
located inside the ergosphere, and our numerical simulations 
are fully consistent with this result (see fig. 4, Al). All these 
allow us to conclude that there is no causality clash associ- 
ated the electrodynamic mechanism in general and with the 
Blandford-Znajek solution in particular. 



7 CONCLUSIONS 

• The 3+1 equations of black hole electrodynamics can be 
written in more general and simpler form than the classical 
formulation of Thorne and Macdonald (Thorne & Macdon- 
ald 1982; Macdonald & Thorne 1982; Thorne et al.1986). 
Our equations are very similar to the classical equations of 
electrodynamics in matter and hold in any system of co- 
ordinates which delivers time-independent metric form, e.g. 
Kerr-Schild coordinates. 

• The inertia of magnetospheric plasma is not an essential 
element in the mechanisms of extraction of rotational en- 
ergy of black holes involving electromagnetic field. We have 
seen no indications of any major deficiency of the electro- 
dynamic model which would suggest a need for the more 
complex MHD model, at least at the black hole end. The 
force-free approximation, however, may break down locally 
via formation of dissipative current sheets. We have con- 
structed a simple model for radiative resistivity based on 
the inverse Compton scattering of background photons and 
used this model to show the development of the equatorial 
current sheet within the ergosphere of a black hole with a 
dipolar magnetic field. This current sheet supplies energy 
and angular momentum for the surrounding force-free mag- 
netosphere. 

Having said that, we need to stress that the role of particle 
inertial has not been considered in this paper and requires 
further investigation. For example, one would expect the 
inertial effects to become important in the remote part of 
the outgoing wind, like they seem to do in the winds from 
neutron stars, e.g. (Mestel 1999; Kennel & Coroniti 1984). 

• The gravitationally induced electric field is the ultimate 
cause of the poloidal currents in the black hole magneto- 
spheres. Within the ergosphere this field cannot be screened 
by any static distribution of electric charge. This makes a 
"magnetized" black hole very different from the classical 
unipolar inductor (or any other battery), where the poten- 
tial difference between terminals is created and sustained by 
the electromotive force which first separates electric charges 
and then drives electric currents against the electric field. 
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In fact, within the ergospheric current sheet the cross- field 
current flows along the electric field but not against it. The 
special role played by the ergosphere allows us to call it the 
"driving force" of the Blandford-Znajek mechanism. Since 
the same applies to the otherwise rather different Penrose 
mechanism (Penrose 1969), this suggests that it is the exis- 
tence of the ergosphere that makes possible energy extrac- 
tion from a black hole in any form. 

• The Blandford-Znajek solution in particular, and the 
electromagnetic mechanism of extraction of rotational en- 
ergy of black holes in general, do not clash with causal- 
ity. First of all, the Znajek horizon condition is not a 
boundary condition but a regularity condition at the fast 
critical point. This perfectly justifies its utilization in the 
Blandford-Znajek monopole solution. Secondly, the ergo- 
sphere is causally connected to the outgoing wind by means 
of both fast and Alfven waves. This conclusion is strongly 
supported by numerical simulations that show that the 
Blandford-Znajek monopole solution is asymptotically sta- 
ble. 

• Our results fully agree with the conclusion of Punsly 
and Coroniti (Punsly & Coroniti 1990a; Punsly & Coroniti 
1990b) and Beskin & Kuznetsova (2000) on the superficial 
nature of the interpretation of the event horizon as a unipo- 
lar inductor of black holes. The failure of the Membrane 
paradigm to predict the outflow of energy and angular mo- 
mentum along all magnetic field lines penetrating the er- 
gosphere has clearly exposed its limitations. Although the 
analogy with conducting sphere is based on mathematically 
sound grounds, it does not account for all important prop- 
erties of black hole electrodynamics and does not reveal the 
true physical nature of the Blandford-Znajek mechanism. 
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APPENDIX A: LIGHT SURFACES 

The physical significance of the light surfaces in the black 
hole magnetospheres is well understood and indisputable. 
However, the author of this paper could not find any more or 
less complete and systematic description of their properties. 
This section is intended to fill the gap and to provide a 
convenient reference for the reader of the paper. 

Al Basic properties of light surfaces 

Consider a point orbiting a Kerr black hole with angular 
velocity SI. The type of its world line depends on the sign of 
the function 

f{Sl,r,6) = g^S? + 2g t4 ,Sl + g tt . (Al) 
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It will be time-like (subluminal rotation) if / < 0, space-like 
(superluminal rotation) if / > 0, and null if / = 0. Similarly, 
a magnetosphere, rigidly rotating with angular velocity SI, 
is divided into the regions of subluminal and superluminal 
rotation depending on the sign of f(Sl,r,0). The separating 
surfaces, defined via 

/(n,r,0)=O, (A2) 

are called "the light surfaces" . Here we derive the basic prop- 
erties of light surfaces of magnetospheres with < SI < Slh, 
where Slh = a/2r+ is the angular velocity of the black hole 
(without any loss of generality we assume that a > 0.) 

(i) The light surfaces are symmetric relative to the equa- 
torial plane. 

This property immediately follows from the symmetries 
of the metric tensor. 

(ii) There exist only two light surfaces, the inner one and 
the outer one. Inside of the inner surface and outside of the 
outer surface the rotation of the magnetosphere is superlu- 
minal. 

Consider the case SI = 0. Then 
f(Sl,r,6) = g tt (A3) 

and, thus, the only light surface is the ergosphere of the 
black hole. Because gtt is positive inside of the ergosphere 
and negative outside of it, the rotation inside of this light 
surface is superluminal. Since f(Sl,r,8) is analytical, this 
light surface will continue to exist for SI > 0, though its 
location may differ. We shell call it the "inner light surface". 
Since along the rotational axis eq.(A2) is reduced to g tt = 
for any SI, this surface always includes the point (r, 6) = 
(r+,0). 

If < fi <C 1 then another light surface emerges from 
infinity. Indeed, for r 1 

f(Sl,r,9) = r 2 sin 2 6S1 2 - 1 (A4) 
and, thus, eq.(A4) has the solution 

r =nkre' (A5) 

which moves to infinity as SI — > 0. We shell call this light 
surface the "outer light surface". From (A4) it follows that 
outside of this surface the magnetosphere rotates superlu- 
minally. Notice, that according to eq.(A4) no other light 
surface can continue to infinity for any value of SI. 

Since f(Sl,r,6) is analytical, an additional light surface 
can only 1) bifurcate from a point where f(Sl, r, 6) has a 
local minimum or maximum, or 2) appear as a result of 
splitting of already existing light surfaces. To exclude the 
first option we notice that 

p 4 df/d8 = [(Ap 4 + 2r(a 2 + r 2 ) 2 )^ 2 - 

-4ar(a 2 + r 2 )Sl + 2a 2 r] sin26>, 

which is positive for r > r+ with exceptions of the rotational 
axis and the equatorial plane where it vanishes. However, 
the only point of the rotational axis where f(Sl, r, 9) ever 
vanishes belongs to the inner light surface. Similarly, there 
are only two locations on the equatorial plane where / ever 
vanishes and those belong to the inner and the outer light 
surfaces for all values of SI. Indeed, along the equator / = 
s(Q, r)jr where 
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s(A, r) = (r 3 + ar + 2a 2 )fl 2 - 4afl + 2 - r. 



(A6) 



Since s is a cubic polynomial in r, it has at most three roots. 
The extremes of the cubic are given by 



(A7) 



which shows that one of the roots is always negative. For 
fl <iC 1 the two other roots are r in — 2 and r out = 1/fl. 
Obviously they belong to the inner and the outer light sur- 
faces. This also proves that no additional light surfaces can 
appear via splitting of the inner or the outer ones. 
Finally, on the event horizon 

f(Q,r,9) = ^± sin 2 9(Q-Q h ) 2 , (A8) 
P+ 

showing that no light surface can bifurcate from the event 
horizon either. 

(iii) Inside the inner light surface the rotation is "too 
slow" whereas outside of the outer light surface it is "too 
fast". 

The condition of subluminal rotation, /(fl, r, 9) < 0, leads 
to the following constraint on fl: 



(A9) 



where 



= n z - 



flz + - g~ttfgH>, 



fln 



At) 



where Q. z = —gt<p/g<j><i> is the angular velocity of the local 
zero angular momentum observer (ZAMO). On the horizon 
fl < Qh — Qmin and by continuity fl < fl m j n in the whole 
of the inner superluminal region. Similarly, for r —> oo one 
has Q, m ax < 1/r sin 9 < fl and by continuity fl > fl max in 
the whole of the outer superluminal region. 

(iv) The inner and the outer light surfaces do never have 
common points. 

Since, as we have established above, df/d9 > in the 
space between the rotational axis and the equatorial plane, 
all we need to show is that the root r in of eq.(A6) is always 
less than the root r out . We already know that this is true for 
fl <JC 1 and, hence, it is sufficient to show that these roots of 
eq.(A6) never merge. From eqs.(A6,A7) one finds the merger 
condition 



£(a,fi) 



(1 +afl) 3 +27afl 3 



27fT = 0. 



It is easy to see that d^/da > 0. Moreover, direct calcula- 
tions show that £(a, fl/^a)) > 0, where the equality holds 
only if a = 1. These ensure £(<z, fl) > for < fl < fi/,(a). 

(v) The inner light surface is located between the ergo- 
sphere and the event horizon. 
At the ergosphere 

f(Q,r,9)=g^a(n-2Q z ). 

For < Q <C 1 one has fl < Q. z and / < 0. Thus, the inner 
light surface shifts inside the ergosphere as ft turns positive. 
A light surface can also make contact with the ergosphere 
when Q — 2Q Z but that can only be the outer light surface. 
Indeed, simple calculations show that dfl z /dr < for r > 
r+. Thus, along the radial ray passing though the point of 




Figure Al. The light surfaces (thin lines) of a rigidly rotating 
magnctosphere. Here Q = 0.7Qh and a = 0.9. The thick line 
shows the ergosphere. 



contact /(fl, r, 9) is positive all the way from the ergosphere 
to infinity. 

The conclusion that the inner light surface always remains 
outside of the event horizon simply follows from eq.(A8). 

(vi) The coordinate r of the inner light surface monoton- 
ically increases with 9 for < O < it/2. For the outer light 
surface there holds exactly the opposite. 

Direct inspection of the lights surfaces found in the limit 
fl <JC 1 (see eq.A3 and eq.A5) immediately shows that they 
have this property. Thus, one need only to verify that dr/d9 
does not vanish anywhere on the critical surfaces with the 
exceptions of the equatorial plane and the symmetry axis. 
Since 

dr^ dl d£ 
d9 09 1 dr 

this follows from the fact that df/89 does not vanish for 
< O < it/2 (see the proof of property (ii)). 

The geometrical properties of light surfaces are summarized 
in figure Al. 



A2 Light surfaces and magnetically driven winds 

Consider a force- free region of a steady-state magnetosphere. 
As we have established above, each magnetic field line that 
penetrates the black hole ergosphere is forced to rotate with 
angular velocity < fl < flh and may cross the light sur- 
faces corresponding to this value of fl. Let us show that in 
the inner superluminal region charged particles are forced 
to move towards the black hole whereas in the outer super- 
luminal region they are forced to move away from it. To 
simplify calculations we shell use the Boyer-Lindquist coor- 
dinates and assume, without any loss of generality, that the 
magnetic field line under consideration is outgoing (B r > 0). 
Any such particle slides freely along the magnetic field 
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line and participates in the drift motion across the field line. 
Thus, its velocity vector is 



v = D xB/B +kB, 



(A10) 



(for simplicity sake we ignore particle's gyration.) The radial 
component of v is limited from above and from below via 
the condition v 2 < 1. In the limit v 2 = 1 one has v = v±, 
where 



v± = ^(DxB±BVB 2 -f 2 ) 



(All) 



Using eqs. (29,75) to express B 2 — D 2 and DxB in terms 
of B and fl one obtains the following result for the poloidal 
component of v. 



v 



p± V± Bp , 
where 



(A12) 



B 2 

73/(^,0) 



(n - 

Given the assumed direction of B p , v p+ has the highest 
possible and v p _ has the lowest possible radial component. 

Since in the interior of the inner light surface Q < £l z 
and / > we conclude that r) + is negative there (notice 
that H,f> < 0) and, thus, charged particles are bound to slide 
towards the black hole. As the result, no magnetic field line 
within the inner light surface can turn away from the black 
hole unless this takes place in the region where the force-free 
approximation breakes down (e.g. within the equatorial cur- 
rent sheet of the magnetospheric Wald problem, see figure 
40 

Similarly, one finds that in the outer superluminal re- 
gion ij- is positive and, thus, charged particles are forced 
to slide away from the hole (Basically the same analysis was 
used by Goldreich & Julian (1969) with application to pulsar 
winds.) 



A3 Light surfaces and Alfven critical surfaces 

The normal wavespeeds of Alfven waves in force-free elec- 
trodynamics in FIDO's frame are given by 

fj,± — ■ n (A13) 

where 

/i ± = (l)xB ±B^fB 2 - D 2 ) IB 2 (A14) 

(Komissarov 2002a). By definition, at a critical Alfven sur- 
face with the outgoing unit normal n either fj, + or p,- 
changes its sign. As pi± are given by exactly the same ex- 
pressions as v± (All) we immediately obtain that 



p± = (n-B p )r)±. 



(A15) 



Thus, in the force- free limit the Alfven surfaces coincide with 
the light surfaces. If B p is outgoing then p+ corresponds to 
the outgoing wave and vanishes at the inner light surface. 
Inside this surface p, + < and, thus, no Alfven wave gener- 
ated in this region can escape from it. Similarly, any Alfven 



wave generated in the exterior of the outer critical surface 
is bound to remain there. 

Koide (2003) claimed that in the force-free limit the 
inner Alfven surface coincides with the event horizon. This 
is obviously incorrect due to the property (v) of Appendix 
Al. 



APPENDIX B: BOYER-LINDQUIST AND 
KERR-SCHILD COORDINATES 

In the Boyer-Lindquist coordinates {t, <f>, r, 8} the Kerr met- 
ric is 



ds 2 = g tt dt 2 + 2g, 



+ 



2 + g rr dr 2 + geedO 2 , (Bl) 



where 

gtt = z - 1, g t <i, = -zasm 2 6, g^ = E sm 2 d/p 2 , 



2 / a 2 

g rr = p /A, goo = p , 



and 

2 

P = 

Z — 

S = 

A = 



r 2 + a 2 cos 2 0, 
2r/p 2 , 

I 2 i 2\2 2 A • 2 a 

(r + a ) — a Asm 0, 

2,2 r, 
r + a — Ir. 



Notice that we use such units where the black hole mass 
M — 1. At the event horizon, where A = 0, one has g rr = co 
which is the well know coordinate singularity of the Boyer- 
Lindquist coordinates. The shift vector is 



P = (St^/5^,0,0), 



(B2) 



so the Boyer Lindquist FIDO is in the state of purely az- 
imuthal motion with constant angular velocity. From (3) 
one finds that = and, thus, FIDO is also ZAMO (Zero 
Angular Momentum Observer, Bardeen at al.,1973.) At the 
even horizon its world line turns space-like and this is an- 
other consequence of the coordinate singularity. This coordi- 
nate singularity can be removed via the following integrable 
coordinate transformation, also singular at the event hori- 
zon: 



dt - 
d<)> - 

where 



dt + G(r)dr, 
d<)> + H(r)dr, 

G(r) 



-2r/A, H(r) = -a/A. 



(Since the transformation law for t depends on r this intro- 
duces a different foliation of space-time.) The corresponding 
transformation matrices are 



(A 



i\p 

!p> 



1 





G 


\ 






1 


H 












1 















1 J 






<\ 





-G 









1 


-H 





■ 







1 















1 



(B3) 



(B4) 
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Here we assume that 1) the primed indexes refer to the 
Boyer-Lindquist coordinates; 2) the upper index indicates 
the matrix row and the lower index indicates the matrix 
column. This allows us to apply the Einstein summation 
rule to the transformation laws, e.g. 

In the new coordinates, known as the Kerr-Schild coor- 
dinates, the Kerr metric is 



(B5) 



ds — gttdt + 2g t ^dtd(f> + 2g tr dtdr+ 

g</,<t,d4> 2 + 2grtf,dcj>dr + g rr dr 2 + geedO 2 , 



where 



g u = z-\, g t <t> = -2a sin 9, g tr = z, 
g^ = Esin 2 9/p 2 , g r4> = -asin 2 6(1 + z), 

g rr = 1 + z, gee = p . 



Thus the spatial metric tensor is 

/ Ssin 2 (6»)/p 2 -asm 2 (6)(l + z) 
Jij = -asm 2 (6)(l + z) 







1 + z 




(B6) 



As g rr f, 7^ the spatial coordinates are no longer orthogo- 
nal. However, at infinity this tensor has exactly the same 
form as the metric tensor of Euclidean space in spherical 
coordinates. This is no surprise as for r — » co one has 
F(r),H(r) — > and the Kerr-Schild coordinates become 
identical to the Boyer-Lindquist coordinates. 

The Kerr-Schild lapse function and its shift vector are 



a = 1/Vl + z, 

z 



(0, 



1 + z 



0). 



(B7) 
(B8) 



Thus, the Kerr-Schild FIDO is moving radially towards the 
true singularity at r = 0. 

It is illuminating to determine the motion of the Kerr- 
Schild FIDO in the Boyer-Lindquist coordinates. Simple cal- 
culations give 

v 4, = 2ar/£, v r = — 2rA/E, v 6 = 0. 

Thus, the Kerr-Schild FIDO has the same angular veloc- 
ity as the Boyer-Lindquist FIDO but also moves radially 
towards the true singularity. 

The tensor transformation laws 



and 



allow us to relate the components of the four-potential and 
the electromagnetic field tensor in both coordinate systems. 
It is easy to see that all components of the four potential 
except U r and all components of F^ v with p, v 7^ r are in- 
variant. This immediately shows that the scalar potential 
$ = — lit-, the magnetic flux function ip = = "Ft^ 

and = F^, t are the same in both coordinate systems. 
Since the angular velocity of magnetic field is given by 

Q = -Fte/F^e, 

it is also invariant. 



APPENDIX C: NUMERICAL METHOD 

CI Augmented system of the method of 
generalized Lagrange multiplier 

From eq.(15) one finds 

d t (V-B) = 0. 



(CI) 



This well known result shows that it is sufficient to enforce 
the divergence free condition (14) only for the initial so- 
lution and it will automatically be satisfied at any t. Un- 
fortunately, straightforward applications of many numerical 
schemes, perfectly suitable for other hyperbolic systems of 
conservation laws, fail to deliver a good result in the case 
of electrodynamics and MHD simply because their discrete 
equations are not consistent with any discrete analogue of 
(CI). In particular, this applies to the method of Godunov 
which has many beneficial properties and is currently consid- 
ered as generally superior to many other numerical schemes 
for hyperbolic systems. There have been many attempts to 
find a cure for this "div-B problem" , though no perfect so- 
lution seems to have been found so far ( see the review in 
Dedner et al.,2002.) 

One of the ways to handle this problem involves con- 
struction of a somewhat different system of differential equa- 
tions, the "augmented system", where the divergence free 
condition (14) is no longer included and V-B may be trans- 
ported and/or dissipated like other dynamical variables. The 
idea is not to enforce the divergence free condition exactly 
but to promote a natural evolution of the system towards 
a divergence free state. Provided the augmented system is 
hyperbolic, it can be solved numerically using the Godunov 
method (Munz et al.2000; Dedner 2002). 

Following this idea we propose to modify eq.(8) as fol- 
lows 



V *F a,3 + V Q * + Kfc Q * = 0, 



(C2) 



where Vl/ is a scalar field which we shell call "the pseudo- 
potential" and k =const. As before k" — dt, where t is the 
global time coordinate of space-time. From equation (C2) 
one immediately obtains the evolution equation for * 



□ * + KV t * = 0, 



(C3) 



where □ = V^V" is the d'Alembert operator. This is the 
well known "telegraph equation" and it prescribes both 
propagation and dissipation of >F In fact, V B is governed 
by a similar equation. Indeed, from the time component of 
eq.(C2) one obtains 

V** + = -V B, 

a 

whereas from the spatial part one has 

v« (^) =-v t v l *. 

Combining, these two equations one finds first 



v*v t (-^— ) + «v t (-^— ) + V t V l (V** + k#) = 

id, finally, 

(™)+-(™)=°> 
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Direct 3+1 splitting ol eq.(C2) gives us two conservation 
laws: 

ft (^D+Mv^' -§/?')) = 

(C5) 



*(«a^-«<(>/7^)) 



(C6) 



and 

ft (V7(B l + + d, {Vl(e ijk E k + a*g ij )) = 

Summarizing, we have constructed a system of two vector 
conservation laws, equations (57, C6), and one scalar conser- 
vation law, equation (C5). 



C2 Numerical scheme 

All 7 evolution equations of the augmented system, as well 
as the equations of the general relativistic electrodynamics 
itself, are conservation laws. We can write them as a single 
abstract vector equation of the form 

dt(^yQ K ) + ft i^ Kj ) = yfrisfa + 



(C7) 



where 



B i 



(C8) 



are the conserved variables, where i is the index of spatial 
coordinates, 



T Ki = (s 3 - e ijk E k + a^g",-e ljk H k ^ 

are the corresponding hyperbolic fluxes, 



(C9) 



S (g) ~ 



-0 



-f9*(>/7— ), -^ft(V7«3 ),p(-a«S + /3 1 ) (CIO) 



\/7 « 
are the geometrical source terms, and 

Sfc = (aK<f,0\^a( n Dl+a ± DY)) (Cll) 

are the potentially stiff dissipative source terms. 

To handle stiff source terms we use the time-step split- 
ting technique (LeVeque 1997). Namely, each time-step we 
first integrate the truncated system 



d t (V^Q K ) + ft(^ Kj ) = 



(C12) 



The result is then considered as the initial solution for an- 
other truncated system 



d t Q K = S K 



(d). 



(C13) 



which is integrated over the same time step. 

In order to integrate (C12) we apply the Godunov 
method, slightly modified to accommodate for the space- 
time curvature (Pons et al.1998). Integrating (C12) over a 
spacetime cell (At x Aa; 1 x Ax 2 x Aa; 3 ) one obtains 

Q K (t + At) = Q K (t) - H JV^'AS, + S&At, (C14) 

i 

where the summation is taken over all spatial faces of the 
cell. In this equation Q is the vector of conserved variables 
averaged over the cell volume, 5( 9 ) is the vector of source 
terms averaged over the cell volume and the time interval, 
and T Ki are averaged over the corresponding j-face of the 
cell and the time interval. If the computational cells are 



created by coordinate surfaces, as it is in our case, then 
AV = v /jAx 1 Ax 2 Ax 3 and AS Z = ±^yAx j Ax k , where 
i^j^k. 

Equation (C14) allows us to advance the numerical so- 
lution in time provided T Kj and S( g ) are known. In the Go- 
dunov method, the fluxes are found via solving the cor- 
responding Riemann problem at the cell interfaces (LeVeque 
1997). Following Pons et al.(1998) we setup such problems 
utilizing the "primitive" solution vector, P, which includes 
the components E* = D l and B % of electric and magnetic 
field as measured in the orthonormal basis of local FIDO, 
P K = (*, B\ D 1 ) 1 . This allows us to consider the Riemann 
problem as observed in the locally inertial frame of the FIDO 
initially placed at the cell interface and use a special rela- 
tivistic Riemann solver to find its solution in this frame. 
At this stage, it is important to take into account the finite 
speed of the interface relative to FIDO , ji % (Pons et al.1998). 

In the locally inertial frame with the Riemann disconti- 
nuity normal to the x-axis the relevant system of ID equa- 
tions is linear 



d t P K + Afd x p L 



o 



(C15) 



with the Jacobean matrix 



A 






















1 



1 



-1 



The eigenvalues of this matrix 

M(i) = M(2) = M(3) = 1; 

M(4) = 0; 



M(B) — M(6) — A*(7) 



= -1 



provide the wavespeeds of the hyperbolic waves. Other prop- 
erties of these waves are given by the eigenvectors of A. The 
right eigenvectors are 

r (1) = (1,1,0,0,0,0,0) 

r (2 ) = (0,0,-1,0,0,0,1) 

r (3) = (0,0,0,1,0,1,0) 

r (4) = (0,0,0,0,1,0,0) (C16) 

r (5) = (-1,1,0,0,0,0,0) 

r (6) = (0,0,1,0,0,0,1) 

r (7) = (0,0,0,-1,0,1,0) 

Since A is symmetric, its left eigenvectors, i& , coincide with 
the corresponding right eigenvectors: 



(0 



K 

' r W> 



(C17) 



It is easy to see that the solutions 2,3,6, and 7 describe 
the usual electromagnetic waves. The solution 4 simply re- 
flects the fact that all waves of vacuum electrodynamics are 
transverse and any discontinuity in the normal component 
of electric field is due to a surface electric charge distribu- 
tion, the solutions 1 and 5 describe new waves which do not 
exist in electrodynamics; they transport Vl' and V • B with 
the speed of light. 
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The solution to the Riemann problem with the left and 
the right states, P( ( ) and P( r ), and the interface speed f3 x 
is 



P = < 



where 



^(0 

P(r) 

'(0 + E K w r w 

»=1,3 



if P x < -1; 

if P x > +i; ( ci8) 

if -1< f3 x < +1, 



t(i) 



(P (r) -P (l) )-r 



(0 



"(0 • r w 



Notice, that the 4th wave is ignored and the x-component 
of D is set to be the mean value of the left and the right 
states 



D t = 0.5(4, +Df r) ). 



(C19) 



In order to make the scheme second order in space and 
time we introduce slope-limited linear distribution of primi- 
tive variables within each cell and use the half-time-step so- 
lution following the description in (Falle 1991; Komissarov 
1999). 



1.3 - 
Bz 
1.2 - 

1.1 - 



-1.0 



0.0 



1.0 



Figure CI. Stationary Alfven wave. The exact solution and the 
initial numerical solution are shown by the continuous line. The 
crosses show the numerical solution at t = 2 



C3 Test simulations 

To test the code we carried out a number of test simulations. 
Here we describe only two of them. In both cases the space- 
time metric is Minkowskian, so E = D and H — B. 

The dissipative time scales of our resistivity model is 
tied to the computational time-step, At. Namely, we set 



d/At, 



(C20) 



and, following the analysis of Sec. 2. 3, we make the cross field 
resistivity strongly dependent on B 2 — D 2 : 



ox 



if 

b(D ± - DD/Dl if 



B 2 > D 2 
B 



2 < jj2 



(C21) 



where (D*l) 2 = B 2 — D 2 . This model provides continuous 
a± at D 2 — B 2 which is preferable for numerical reasons. 
The corresponding solutions to equations (C13) are 

-CTllt 



£>ll(t)=I>||(0)e- 



and 

D±(t) 



D1 + 



D1(D±(0) - Dl)e- b °l 
Dx(0) + (£>!- £>x(0))e" 



(C22) 



(C23) 



In fact, each time-step we first evolve By as in eq.(C22) and 
then D± according to eq.(C23). In most of the problems we 
use d < 1 and b = 0.1. In addition, we assume that the drift 
current is described by eq.(100). Although this model resis- 
tivity is too high compared to the physical one, any much 
lower model resistivity would not allow to resolve current 
sheets (and would become lower that the numerical resistiv- 
ity.) 



C3.1 Alfven wave 

The system of force-free electrodynamics allows two types 
of hyperbolic waves: 1) the fast waves, which have the same 



properties as linearly polarized waves of vacuum electrody- 
namics, and 2) the Alfven waves, which are similar to the 
Alfven waves of relativistic MHD (Komissarov 2002a). The 
analytical solution for the Alfven wave used in this test de- 
scribes a stationary wave with B x — B y — D z = 1, D y — 0, 



B z (x) 
and 



1 + 0.15(1 + sin57r(a; - 



1 

0.1)) 
1.3 



for 
for 
for 



x < 0; 

< x < 0.2; 

x > 0.2; 



D x — —B z . 



The computational grid included 200 identical cells equally 
spaced between x — —1.5 and x = 1.5. 

Figure CI shows the solution at t = 1. While an electro- 
magnetic wave would cover a distance Ax — 1 during this 
time, the Alfven wave remains stationary and only spreads 
due to finite resistivity. 

C3.2 Current sheet 

Here we consider the following Riemann problem: 
D = 0, B z = 0, B x = 1 



Bo 
-Bo 



for 
for 



x < 0; 
x > 0; 



The symmetry of this problem implies B y — at the in- 
terface at time t > 0. If Bo < 1 then this problem has a 
force- free solution describing two fast waves, switching-off 
the tangential component of magnetic field and switching- 
on the z-component of electric field (up to the value of —Bo.) 
In fact, the vacuum electrodynamics has exactly the same 
solution to this problem. The numerical solution at t = 1 
for Bo = 0.5 is shown in the panel a) of figure C2. The 
computational grid is uniform with 200 cells. 
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Figure C2. Riemann problem described in Section C3.2. Left panel: Numerical solution at t = 1 for Bo = 0.5. Right panel: Numerical 
solution at t = 1 for Bq = 2.0. The current sheet is located at x = 0. 



If Bo > 1 then the vacuum electrodynamics still allows 
solutions of this kind. However, the force-free electrodynam- 
ics has no solution in this case because B — D vanishes 
when \B y \ drops down to the value of (B + l)/2Bo- At this 
point the fast wave is terminated by increased cross-field 
conductivity that locks B 2 — D 2 to 0. Downstream of these 
waves, the Poynting flux is directed towards the interface lo- 
cated at x = where the electromagnetic energy dissipates 
in a current sheet. The panel b) of figure C2 shows the nu- 
merical solution for Bo = 2 at time t = 1 on a uniform grid 
with 200 cells. As expected, the global solution is not sensi- 
tive to the exact value of the cross field conductivity which 
only effect the structure of the current sheet, the higher t± 
leading to the wider current sheet. 
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